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ABSTRACT 


The advent of the missile and satellite era has brought 
with it the need for large and highly accurate radar tracking 
systems. Thus, the antenna tracking error generated by wind 
disturbances has become a significant consideration. This 
study presents a compensation technique which substantially 
reduces these disturbance errors without degradation of the 
system response to the tracking input. 

The study was proposed by Philco Corporation, Western 
Development Laboratories, Palo Alto, California. The proposal 
resulted from tests conducted on a production model of a 
satellite tracking system having a 60 foot, solid skin, high 
inertia antenna. A large portion of this thesis was accom- 
plished at the Philco Laboratories during the Summer of 1960. 

The authors wish to express their appreciation to 
Dr. M. P. Pastel of the U. S. Naval Postgraduate School for 
his assistance and counsel in the preparation of this thesis 
and to Mr. L. J. Harvey of Philco Corporation for his assistance 
in the development of the problem and for the generous facili- 


ties made available to the authors. 
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i limsroduction 

The purpose of this study is to present a compensation 
technique for reducing antenna tracking errors caused by 
wind-gust disturbances. The antenna system upon which this 
study is based is the position control system for a 60 foot 
radar antenna. This system meets. the desired response to the 
tracking input. However, the tracking error caused by wind 
gusting is not acceptable. The compensation utilized must 
reduce this error without affecting the system response to 
the tracking input. This problem corresponds to having an 
antenna system installed at the tracking site which meets 
the designed input tracking response. Data showed excessive 
tracking error resulting from dynamic wind pressure on the 
solid skin of the large, high inertia dish. This error must 
be reduced by adding some type of compensation that will not 
modify the designed tracking input response. 

An error of a few seconds of arc can be excessive in 
determining the precise position of a high altitude orbiting 
satellite. To meet such stringent error regulating require- 
ments while, at the same time, achieving acceptable tracking 
input response using conventional compensation techniques 
poses a formidable problem. Compensating techniques require 
the sensing of the effects of the wind loading on the antenna 
system. Wind sensing might take the form of external devices 
such aS sensitive annamometers. These might be placed in 
locations around the periphery of the antenna and at distances 


so as to determine in advance the load disturbance which the 





antenna would see. Although the nature of a prediction 
sensing device lends itself to the development of the ideal 
regulation system, the accuracy of the prediction must of 
necessity be at least as accurate as the tolerance of the 
desired output response. To date, no wind sensing trans- 
queers have been Peeduced that measure accurately transient 
wind conditions over a sufficiently ean aiiérement of time. 
Even if such devices were feasible and available, their 
sensing would not necessarily reflect the actual loading 

of the antenna dish some distance away. A mountainous site 
might be expected to show large drafting and eddy-currents 
about the dish which would preclude the accuracy of any 
external sensing device to determine actual dish loading. 

A second and more reasonable approach might take the form 
of a sampling of the torque on the antenna drive shaft. 
Corrections would then be made in the control system based 
upon these samplings. This approach is complicated and 
expensive, and at least part time use of a digital computer 
would be required. 

A simpler and more direct approach would be to continuously 
compare the output position under external loading with the 
output obtained having no external load with the same tracking 
input. The difference of these outputs would then be a 
measure of the effects of transient wind loading on the 
antenna. Although this approach in no way predicts the loading 
effect, it does permit the accurate measurement of the effect 
which in turn may be compensated at a faster rate than that 


possible by the basic system. 





The solution presented here will be developed in the 
above manner using as the basic approach the conditional 
feedback concept set forth in the American Institute of 
Electrical Engineers' publication, "Applications and 
imow@stry,'' July 19e5 > by G. Lang and J. M. Ham. The condi- 
tional feedback system permits design specifications on 
input-output response and on disturbance-output response to 
be met independently. 

To prevent repetition and to simplify the presentation, 
only the azimuth axis of rotation will be considered. Design 
criteria for the compensation will be developed based on 
selected arbitrary requirements for disturbance-output 
response improvement. A specific compensation will be designed 
using the developed criteria and the response will be 
determined for various disturbance inputs. A design technique 
is developed which facilitates the determination of this 
compensation. 

An investigation of the characteristics of a severe wind 
gust disturbance will be utilized to develop a disturbance 
input which will subject the compensated system to the most 
demanding conditions expected. 

An analog computer simulation will be made of the complete 
compensated system. Data will be displayed showing the 
effects of the selected compensation for various wind 


disturbance inputs. 





2. Development of the basic approach. 

In the conventional single-loop feedback control system 
the input-output response, as well as the disturbance-output 
response, are dependent on the feedback. Thus, in most cases 
utilizing conventional feedback control, a compromise of the 
input response must be accepted if a satisfactory disturbance 


response is to be achieved. Fig. 2-1 shows such a conventional 


Geucrol system 





Figure 2-1. Diagram of conventional control system. 


icaerig. 2a G, is the Laplace transform of some control system; 


G., represents the feedback compensation. R is the normal 


Z 
input signal to the system; U is the disturbance signal at 
the output of the system; C is the final system output. 

ie transfemmaunctions associated with Fig. 2-1 are readily 


developed as follows: 


C = G,CR - GC) + 


CCl + G,G,) = GR + U 
Gy 1 
C= ———— R+ —— (2-1) 
1 + G1G>5 1 + GG, 


Thus, the feedback G, affects both input and disturbance 





responses. However, the problem treated in this thesis is 
one in which the effects of disturbance must be compensated 
without compromise of the given input-output response. A 
compensation technique meeting these requirements has been 
presented in "Conditional Feedback Systems" by G. Lang and 
J. M. Ham (1) and will be used as the basis for the approach 
of this thesis problem. 

A basic configuration for a linear conditional feedback 
system is shown in Fig. 2-2, where G,, G., R, U, C, are defined 


as in Fig. 2-1. 


U 
+ 
R 
—>(X) Gy - >} - 
£ 
e 
- 45 
x) H 


Figure 2-2. Basic conditional feedback control system. 


ive output @mer Fig. 2-2 is 


G,(1 + BG.) it 
C= —————-—- 1 R + ———— _U. (2-2) 
1 + G,G5H 1 + G,G.,H 


Let B in equation (2-2) be defined as 


B= GH. (2-3) 


Equation (2-2) then becomes 


C =G)R + ——+~—_—_—_U. (224) 





Petting, U and R be@aitemnately zero, 


(A) 9 mc, (225) 


and 


G | E i : (2-6) 
U R=0 1 + G,G,H 


Equation (2-5) shows that, when B is defined by equation 

(2-3), the input response is dependent only on G,. Thus, G.H 
can be utilized to modify the disturbance response without 
aifecting the input response. The output sensing device 

would normally be described by H. If it is assumed that 

H = 1, then B becomes equal to G, (equation (2-3)). The 
purpose of the function B is shown by considering the system 
operation in Fig. 2-2. The feedback signal (b) is equal to (a) 


since H= 1. With B defined as equal to G (d) is equai to 


1? 
(a) with U equal to zero. Thus, the difference signal (e) is 
equal to zero and no compensation signal is developed. 
Clearly, (d) is equivalent to the desired output response and 
thus B is called a "reference model."' The reference model 

is an undisturbed representation of the transfer characteris- 
tics of the antenna drive system. This conditional feedback 


approach permits the use of the feedback function G. to 


Z 


accomplish the desired disturbance compensation without 


Change to the ommeginal input response: 


3. Development of basic system. 

The basic system investigated in this thesis represents 
the azimuth axis position control system of an actual land- 
based tracking radar having a sixty foot diameter antenna. 

The system consists of a hydraulic drive with tachometer 
feedback and suitable reduction gearing, together with cascaded 
electronic compensation and position feedback loop. 

Prior to any possible investigation of the basic system, 
the hydraulic drive system transfer functions had to be derived. 
This hydraulic system consisted of a servo control value and 
amplifier, a D. C. motor-driven constant displacement pump and 
associated hydraulic drive motor, and reduction gearing. The 
hydraulic system was enclosed with a tachometer feedback loop. 


The block diagram of this hydraulic drive system is shown 







eee. Sai ee a a ee 
: i fs 
S) Servo control Hydraulic fae ik y 
Compensated (x) ® | valve motor N 
error signal & | 
amplifier l 
| G | 
i a a. ee 
Tachometer (GF) 
Figure 3-1. Hydraulic drive system. 


The transfer function (G,) is defined as the joint function 
of the servo amplifier, servo control valve, and hydraulic 
drive motor. Neglecting the tachometer feedback (Ge), the 
output position angle to the reduction gearing is (0,), where 
0, = GX (3-1) 


1 h 
and X is the compensated error signal of the basic system. 
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The block diagram of the hydraulic drive system with normal 
antenna loading and external load torque oe due to wind is 


shown in Pig. 3-2. 





Figure 3-2. Hydraulic system with synthesized loading. 


(3-2) 


where as is the torque due to hydraulic motor inertia (J) 


and is defined as 


_ 2 
i sO, (3-3) 


The load torque (T,), as shown in Fig. 3-2, is 


_ 2 a 
T, = J,s°@, + £, sO, + i = K, (0, Z Q,) (3-4) 
91 = 84 ; (3-5) 
N 


where N is the gear reduction ratio of the hydraulic system. 
Substituting equation (3-5) into (3-4) gives 


_ = 2 
a 1 - Q,) = (Js + f,s)0, eT. (3-6) 
N 





or (3-7) 


K 
ee ees LT 


L 2 iE 2 < 
NCJ, s + fis ¥ Ky) (J,s + fis + Ky) 


Substituting equation (3-7) into (3-6) gives 


Kr Ky 
2, > \ 3 Se 5 
TO NCJ; s + fis + K, ) (J,s + fis + K, ) 


(3-8) 
The flow of oil out of the control valve (Q.) is considered 

to be divided into three parts, that lost to system leakage 
(Q,); that lost because of compressibility of the oil (Qa), 
and that actually used to drive the hydraulic motor (Q. > 

ae 8 LT (3-9) 
The flow out of the control valve is a function of the input 
error voltage (X) within some constant (K.) and is defined as 

une X . (3-10) 
The constant (K.) is defined to include the control valve 
constant as well as the servo amplifier gain constant. Brown 
and Campbell in "Principals of Servomechanisms" [2] define 
the separate oil flows of the system as follows: 

Cream! (3-11) 
where L is the leakage coefficient and P is the pressure 


differential of the system. 


Q. =a eneer, (3-12) 


where V is the volume of oil under compression and B is the 


bulk modulus of the oil. 





Qn = am s®, ; ots) 


where dm is the motor displacement and oO, is shown in 


Pig. 3-2. Dimensions of all constants are listed in the table 


of constants on page 12. 


Making substitutions of these flows into equation (3-9) gives 


= V 
KX = Yam sO, Ge 5S) P. Cs 


The pressure differential (P) of the system is related to the 


torque output of the motor (T,) as follows: 


Ty = dmP. (3-15) 


Substitution of equations (3-15), (3-3), (3-4), into equation 


(3-2) gives 


K -K 
T, = dm P= J 82, + = ees + 7 Q, 
N (J, s + fis + K, ) 


K 
, ee 2 


5 (3-163 
NCJ, s + fis : K, ) 


e 9 
Or 

J as K Ca 5? + fs) 

L ie L 

dm dm Ne 
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2 1 
(J, s ~ f,s * K, ) 


AT 
+ ——_—__—__ nn T. . (3-17) 
dm N (Jy s + fis + Ky ) 


Substitution of equation (3-14) gives 
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2 2 
Jags fomecueys, + £.s) © 
KX = dm SQ, Peeler oy S ) = ‘a +" 5 s : 
ss B dm dm N (J, s 





+ f,s + Ky) 


Ge & os) kK 
B LE 
i: 2 “2 
omienedes + fs 4K ) eae 
If the external torque cr is assumed equal to zero, then, 


solving for ~ in equation (3-18) gives 
1 


2 2 
cs K (J. s ete 
i Vv dm dm N CJ, s +f, s+K, ) 








_L meee. LL La 
amen’ (Ciee’ + £.s + kee 
e L L 


(3-19) 
Since the system stiffness factor (K, ) is relatively large, 


equation (3-19) reduces to 


4 BK, dm / V(J_ + Ji N2) 


. 
Je cp - f, eee, . ee 
Te ox2 a... CCU 
wn Ee? 
N2 N2 Ne 
(3-20) 


The friction coefficient of the antenna drive is considered 


very small with respect to the load inertia constant, thus, 


11 








and equation (3-20) reduces to 


2 2 
N°K_ Boge, + NJ ) 
Vv Ie. m (3-21) 





S s° + LBs + No dm? B 
2 
VJ, + NTI) 
The hydraulic constants associated with this system 


are as follows: 


N = 2160 
K.- game so 5 
V 
Jr = 4.2x10° iboqaeee (referred to antenna) 
Ji = 0.09 ibd noses (referred to motor) 
B = 0.27x10° psi 
E a eo 
DeeeGo. 5x10 ~ in’/Sec/psi 
V =a: 9 aa 
dimmed). 215 in?/radian 


Substituting these constants into equation (3-21) gives 


4 
ce 345x110 (3-22) 


n =a + 42s + 3228) 


Equation (3-22) is the transfer function of the 
uncompensated hydraulic system and was used in the initial 
development of the complete basic system which follows. The 
transfer functions of the complete basic system chosen initially 
for investigation are shown in Fig. 3-3. The values shown for 
the individual transfer functions were suggested by the 
manufacturer aS representing the azimuth control system for 
the antenna under investigation. The open-loop transfer function 


Jw. and the closed-loop function Cr) for the system of 


i 





24 ,050(s+1.7)(st+20) 135,450 


pene eee 
s(st+10) 


s(s2+42s+ 3228) 


(G76086).S 





Go = Cascade compensation 

G, = Hydraulic drive 

G, = Tachometer feedback (tachometer constant = .086) 
N = Gear ratio (2160) 


Figure 3-3 


Block diagram of initial basic system. 


is 





Fig. 3-3 are as follows: 

: 4 

_ oe Cs + 1.7)(s + 20) (3-23) 
: se + HONCS- + 425 + 4038) 


14.96 + oe fier) (Ss +. 20) 


Z 6 6 - 


$ + 3.24x10°s + 5.08x10 


S290 5068 oes + 1.8988x10-s 


(32245) 
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4. Analysis of basic system. 

In order to ascemtain the stability of the System of 
Pee 3-3, a ROWE Seer1terion test was made of the closed- 
loop transfer function OE 


Routh's criterion test: 


1 4.458x10° 3.24x10° 
52 1.8988x10> 5.08x10° 
810 3.14x10° 0 
=1. 1729000 5 .08x10° 0 
3.49x10° 0 0 
5.08x10° 0 0 


The system as given was believed stable, but Routh's 
Criterion test did not verify this assumption. Accordingly, 
a Nyquist plot was drawn of the system transfer function 
in order to observe the degree of instability and is shown 
in Fig. 4-1. Within the accuracy of the plot, the curve 
appeared to pass through the Nyquist point indicating marginal 
stability. For further investigation of the system stability, 
a Bode plot, using asymptotes only, of ae was constructed and 
1s shown in Fig. 4-2. Although the plot crossed the zero db 
line with a negative six db slope, the complex poles of Poo 
having a natural frequency Cw) of 63.5 rad/sec, caused a hump 
in the plot which, if the magnitude were sufficiently large, 
would recross the zero db line and result in system instability. 
In order to estimate the maximum amplitude of the hump, 


that part of the system transfer function (F) resulting from 


i 
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the complex poles can be written in vector form me 


EB a _ ee = Me JN ; (4245) 


e 
Ww : WwW 
(1 = an dS ean 


where: w, = break frequency of the complex poles 
S = damping Coemaeeient 
M = |B 


Nn 
N = phase angle of F 


The maximum absolute value of F ease? is the maximum 

amplitude of the hump. Letting “ = u and substituting into 
Wo 

equation (4-1) gives 


A nt 
— 2. .. 
C1 - u)D+ j2gu 


or 
M = po : (4-2) 


Ci u’) - 43-u- 


Differentiating equation (4-2) with respect to u gives 


dM = =e u’)(-2u) + 8x -u 
du 


7 e 
2 rel = i) + as -u-| 5 fen = u’) + 48 ene | 
(4-3) 


Setting equation (4-3) equal to zero and solving for u gives 


i 
2 


yon 2 -) (4-4) 


Substituting equation (4-4) into (4-2) yields 


‘chestnut and Mayer, Servomechanisms and Regulating 
Systems Design, Vol I, pp 312, September 1953. 
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w 1 
2% (1 - 8 


(4-5) 


max yi 


2 


Equation (4-5) is an approximation of Maes since it was 
developed on the basis that the frequency at Maax was equal 
to the break frequency Cw). However, these two frequencies 
differ slightly. 

The damping factor (8) can be computed from equation 
(3-23) as 


28wo = 42 , 


pmere w, = (4038)2 , or < = ste 


Therefore, from equation (4-5) 


M com. / = 4.62 dp 
max 


The maximum amplitude of the hump, 4.62 db, was added to the 
amplitude at the complex pole break frequency in Fig. 4-2 and 
resulted in a system amplitude at this frequency of 0.12 db. 

Although the above calculation is not exact, it at least 
shows the hump quite close to the zero db line and thus, the 
system has a very small stability margin at best. A decrease 
in system gain would improve the stability by lowering the 
hump below the zero db line, but as only 5.5 db can be dropped 
before the slope crossing the zero db line decreases to -12 db, 
this approach is not satisfactory. 

The above marginal stability was a result of the location 
of the complex poles of Po: These poles in turn were placed 
by the location of the complex poles of the hydraulic drive 


transter Lune tion (G)) and the tachometer feedback (Ge). Near 
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cancellation of these latter poles was achieved by insertion 
of a complex-zero compensator (G.) in series with G,, 
mileviating the G@netepility problem. The design of the complex- 
zero compensator appears in Appendix I. 

Equation (3-22) has shown the complex poles of G,, to be 
ce + 42s + 3228). The form of the transfer function of the 
complex-zero compensator (G.) is shown by equation (I-9) in 
Pependix 1. Thus, the Go necessary to cancel the complex poles 
of G. is 

1 


G, = 3598 SS 


2 4 425 + 3228). 


For the purpose of this study, complete cancellation of the 
complex poles of G, was assumed. In the actual application, 
some residue would remain since equation (1-9) does not 
exactly represent the transfer function of the compensator 
mierie. I-lvor the Appendix. 


Placing Go in series with G, inside the tachometer 


h 


feedback loop results in a new open-loop transfer function, 


myo? (s + 1.7)(sieew@ 


(4-6) 
2 5° (s + 10) 


Fig. 4-3 shows the Bode diagram of the function of 
equation (4-6). A comparison of this Bode with Fig. 4-2 shows 
that the desired cancellation has been achieved without change 
of the system gain or location of the remaining break points. 
Thus, the system response remains the same with the exception 
of the undesirable effect of the complex poles. 


At this point, a limitation was placed on the system 
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cascade compensation (G.). In order that components of less 
stringent design might be utilized in Gos the cascade gain of 
2.405x10° was limited to a value of ioe In order that the 
system gain of equation (4-6) be retained, servo amplifier 
gain of the hydraulic system as well as tachometer feedback 
gain were modified through the following development. 

The open-loop transfer function of the basic system of 
Fig. 4-4 is 


F = ( 


ae (4-7) 


Equating the system gains of equation (4-7) with the desired 


system gain of equation (4-6) gives 


eee) Lae (4-8) 


where K represents the gain constants of the subscripted 
transfer functions of equation (4-7). Substituting oS 1000 
in equation (4-8) gives 

K_K 

; : = o( x2 80 . (4-9) 


h 
KOK LK ¢ 1000 


Substituting Ke = 0.7x0.086, where 0.086 is the tachometer 


constant, into equation (4-9) gives 


(1 - -7x.086)K KK) = 80 


or 


eee o5 ~«C C4=10)) 


ae 





Since all transfer function gains must be positive, equation 
(4-10) indicates that the tachometer gain must be changed in 
order to realize the desired overall system gain. In order 


that K Ky be a positive number, equation (4-10) demands 
(1 - 80K-) 20. (4-11) 


Since the tachometer constant (.086) can not be changed, the 
tachometer gain must be changed so that its value is less 


that 0.1455 in order that KOK be positive. To ensure this 


h 


positive value, a tachometer gain of 0.135 was selected. 


Ke =x .086 ’ 


and 
C1 - -135x.086x80)K K, = 80 , (4-12) 
or 
KOKn = 1150 
is a, 
where K, = 4558 3 therefore 


K. = 1150x3228 = 3.72x10° .. 


h (4-13) 


Equation (3-21) defined K 


N° B am K 


Ne Wie) 1 eee. : 
h V a + N2Im Vv (4-14) 


hoe 


Substitution of equation (4-13) into (4-14) gives 
K = 245 


V 


where K| was previously defined as including the control 


aS 








valve constant as well as the variable servo amplifier gain. 
The transfer functions of the basic system reflecting 


the above changes are 


3.72x10° 

a 2 
s(s~ + 42s + 3228) 
¥ 1 Z 
Ge = 3998 Cs + 42s + 3228) 
Ge = (.135x.086)s 
«ISSO Cs + 1.75 +20) 
c s(s + 10 


The new open-loop transfer function oF) was shown in 
equation (4-6). The new closed-loop transfer function was 


derived: 


_ Je oS + 20> ; CASS > 
© So (s + 1.74)(s% + 45.268 + 724.4) 


This is the basic system transfer function referred to here- 
after as G, and is shown in Fig. 4-4. 

As was shown in the development of the basic approach, 
the conditional feedback compensation effects only the 
disturbance response of the basic system. Pig. 4-5 shows 
the response of the basic system to a unit step disturbance 


input (U), where, 


The response curve of Fig. 4-5 was plotted using the transient 
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response equation determined by the graphical residue 
technique. The details of this procedure are included in 
Appendix II. The maximum overshoot in Fig. 4-5 is 10.3% 
and occurs at t = 0.10 seconds. Fali-time, defined as the 
time interval between ten percent and 90 percent decay during 
the initial cycle of response, is 0.037 seconds. The condi- 
tional feedback compensation to be attempted will try to reduce 
Significantly this fall time. 

Pig. 4-6 shows the frequency response of the basic 
system to a disturbance input of from 0.5 cps. to 40 cps. 
It is noted that as the frequency of the disturbance input 
Signal is increased, the regulating effect of the basic 
system deteriorates. In order to compare this frequency 
response with that resulting after conditional feedback 
compensation, the system bandwidth is defined as that 
frequency above which the input disturbance is attinuated 
by less than 10.62 db, or where M is equal to (1 - .707). 
In Fig. 4-6 this value of amplitude is reached at a frequency 
of 2.10 cycles per second. The corresponding phase angle 


is (-)78 degrees. 
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5. Development of compensation technique. 
The following is a derivation of the conditional 

compensated disturbance-output transfer function 
together with a presentation of a compensation technique. 
Fig. 5-1 is a block diagram showing the inputs and transfer 
blocks of the proposed conditional feedback compensated 
Eyotem, where 

I 


¢. - == the closed loop transfer function of the 


’ 
i ao basic system 


Tl 
I 


the open loop transfer function of the basic system 
= compensator to be determined 
= tracking input information 


= load disturbance information 


G 
R 
7 eS 
C 


= system output 





Figure 5-1. Block diagram of proposed conditional 
feedback system. 
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The transfer function of the output (C) with respect 


to R and U is developed from Fig. 5-1. The relations below 


follow from Fig. 5-1. 


C=U+c=U+#+ (a - C)F. 
c = DF. = (a - C)P. 
b=a-C 

d= 26 c 

e = RG. ae CG. 


a =. Kee 


e=R+ RGoG, =~ CG. 


Combining these relationships gives 


je = se 


(R + RGG, - CG. - C)F. : C5-)) 


Rearrangement of equation (5-1) gives 


(1 + G,G,)F. 


7 U 
C= 7 aE Ec. * | sce (5-2) 
Oo Ores O © 
o 
Samce 1 + EA == > then equation (5-2) becomes 
ah 
7 U ci + GyGo)F, 
— FG — + FG 
Gy 02 G, OZ 
r 1 (l 4 GiGsore 
C= F_ + ? R 
Gg (1 + 6,6.) gl + G65) 
il I 
C = U = eur ($23) 
Girone ol + G_.G_) i a 
O ae 
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Therefore, the output due to load disturbance (U) is: 


i il 
[rae | Lovee, |» “ 


The output due to the tracking input (R) is: 


GR. G5) 


The input-output transfer function Gy is the same as that 

for the uncompensated basic system. Thus, this input-output 

response remains unchanged with the addition of the compensation. 
Root-locus techniques are to be uSed to select a G, which 

will give a desired disturbance-output response. The disturbance- 

output transfer function may be written in the following form 


to facilitate the application of the root-locus technique. 


° 


1 ' : 1 -|; 1 _ non? 
1 + i eet G,G 1 + By + PG DoP> + NONo + ND> 
(5-6) 
where: 
5 
ne 
Oo 
R N 
CG O Oo 
1 1+ RE N + D 
O O Oo 
i "2 
2 D. 
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The denominator of the above function requires two root loci 

in order to determine the final closed-loop roots of the 
disturbance-output transfer function. Since it is desired 

to investigate the effects of several values of the compensator 


G it seemed advantageous if a function, closely representing 


OR 
the disturbance-output function, could be utilized which 
would require only one root-locus plot to place the roots of 
the closed-loop function. An examination of the function 

1 


1 oer G,G, 


shows it to be such a function. This function can 


be rewritten as 


as = 1 + Fo Did, 
1 + GG 1 + Be + B G5 DD, + NN, 


(5-7) 


ee 


where G, = — 


al 
The roots of the function of equation (5-7), hereafter referred 


ney 


to as the "representing function,'' are the same as that for 
the desired disturbance-output function of equation (5-6). 
Furthermore, this function is of a form which requires the 

use of only one root locus in order to place the roots of the 
desired closed-loop function. Thus, if a direct correlation 
between the representing function and the desired function can 
be shown, the problem of selecting G. through simple root- 


locus techniques is made relatively easy. The development 


of this correlation follows. 


SZ 








The representing function of equation (5-7) can also De 


written as 


For the basic type 2 system given, D. takes the form: 


oe 107s + 1) 


1 2 


and represents the denominator of the open-loop transfer 
function of the basic system. 


Assuming the compensator (G,) takes the form: 


m 
Ss (7 4s + 1)(F 5s + 1) 


Cae 38 + Cy Aa 1) 


then the numerator of the representing function, (D.D, + 
N,P>) is some polynomial in s of the form 


(Neos + Bgn-t + . , . eee 


The denominator of the representing function is also of the 
Same general form, (AtsP 4+ Bigee aoe et D'). Since 
no factored s term exists in either the numerator or denominator, 
the response to a step input of this representing function will 
oayeclude a résidwe at the “input pole anserted at the oriein 
and, therefore, will have a steady-state error constant. 
The desired disturbance-output function has been shown 
in equation (5-6) to be: 
D Dd. 


SoS Ne , which is of the form: 
o 2 o 2 o 2 


BS 








s“(+ ahd 1)(7 45s i) (T 38 + 1)(F 4s + 1) 





AtsP + giP-t + . , . . + OD! 


Since this form does not produce a residue at the origin 
of the s plane due to a step input, the time response has 
mo error conS#tanme. 

The time response of the representing function to a 


step input is of the form: 


Ci + Coe a sin(w t +¥), to a second order 
approximation. As time increases, the value of the output 
of the representing function will tend to be larger by an 
amount Cy than that expected for the desired function in which 


wee is equal te zero. 


1 
If there exists no phase difference (0) between the 
representing function and the desired function, then to a 


second order approximation, the correlation would simply be: 


1 1 il 
See u(t) - CO. = | lea Aeuct )*. 
1 +t G,G, 1 1 + ie 1 + GG, 


This, of course, is not the case. Some phase difference 
between the two functions must exist since they are not 
identical. It is this phase difference (% ) between the 
functions which in addition to the constant error term, 
decides the similarity between the two functions. Some 
correlation with respect to this phase difference is necessary. 


An inspection of the desired transfer function, 


oo —————— 
oot a + G,S, DD. + NGP> y NUN> 
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shows that the zeros of this function are made up of the 
poles of FP and G,. An inspection of the representing function, 


( J i Doe i Do? 
1 + G,G9 DOeane Bon 2 . None 


shows this function to have zeros that are made up of the 
poles of G, and G,. Both functions contain zeros which are 
the poles of G,- In addition, both functions have the same 
roots. The only difference between the time response of these 
functions must be due to the difference between the placement 


of the poles of ae and the poles of G (No + Do): the 


1? 
remaining zeros of the two functions. To a second order 
approximation, the phase difference (¥ ) can be determined by 
noting the placement of the poles of Ee and G, with respect 
to the complex pole of either the representing or desired 


function, since they are identical. The phase difference (7), 


determined using graphical residue methods, is as follows: 


% = Sum of the angles from the poles of F_ to 
the complex pole of the representing 
function minus the sum of the angles from 
the poles of G, to the complex pole of 


the representing function. 
(5-8) 


Since the response of the representing function to a 
unit step input, to a second order approximation, is of 
the form 
~X1t 


Cs Coe 


1 sin(w.t +y) , (5-9) 


then the response of the desired function to a unit step 
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input must be, 


ca “i *sin(w.t yer). (5-10) 


3e 
The development thus far permits the determination of 
the correlating parameters, phase difference (%) and steady- 
state error constant (C,)- Techniques for utilizing these 
parameters to predict the desired response from the representing 
response are presented in the following development. 
A comparison of equations (5-9) and (5-10) show that 
they differ by a constant, a proportionality factor, and a 


phase difference (%). The value of Cy can be determined to a 


second order approximation and is shown below. 


UG £2 4? i 
MU vey 
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> =a pole at 
origin 
[-] Poles of R.E. inserted 
for step 
Zeros of R.F. input 
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qe 


Pigure 5-2. Representing function pole and zero placement. 


The representing function, referred to hereafter in this 
section as R.F., has pole and zero locations as shown in 
Fig. 5-2. Applying graphical residue techniques to Fig. 5-2 


for a unit step input gives 
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an Waesisameesto zeros of R.F. from the origin 
divided by the distance to the complex poles 
Of Ree. from the origin. 


Since G. is chosen such that R.F. approximates a second order 


za 
system, the zeros of R.F. are equal to the poles of G,, thus 


he 2 
Ce el 
1 aeCR.F.) 


where W664) is the distance from the complex poles of Gi to 

the origin and WA CR-F.) is the distance from the complex poles 
of R.F. to the origin. Equation (5-11) shows, for example, that 
if the compensated response is to be three times as fast as 


that of Gi, then WA CR-F.) should be approximately three times 


that of W(G,)-. In this example, 


~ (ly? ez 
Ci = (=) ph: 


Thus, for a significant improvement in response time, Cy will 


be small. If the phase difference (3) is also small, C, is 
approximately equal to C, Since both equations (5-9) and 

(5-10) must equal unity at zero-plus time. With both Ci and ¥ 
small the response to a unit step input of the desired function 
will be similar to that of the R.F. For small C, and ¥ the 
overshoot of the desired function may be predicted. The 
overshoot of the R.F. to a second order approximation can be 
determined from the root locus using simple second order 
techniques. The value of C, is then added to this estimate 


of R.F. overshoot and is used as the prediction of desired 


function overshoot. 


oi 





Equation (5-10) shows the phase angle of the desired 
function to be equal to the sum of Y and y . The value of xr 
together with its sign was computed using equation (5-8). 

If a prediction of the desired response based on the R. F. 
response is to be made, some knowledge of Y is necessary 

in order that the effect of * , the phase difference, may be 
known. For example, let ¥ be a Second quadrant angle and ¥& be 
positive. During the initial period of response, the sine 
term of equation (5-10) will take on smaller values than 

will the sine term of equation (5-9). This is to say that 


for small C the response of the actual desired function 


1? 
will be somewhat faster than that of the R. FP. Similar 

examples can be given for negative w or for y in different 
quadrants. Suffice it to say that a knowledge of the quadrant 
location of vy is necessary. The shaded area of Fig. 5-2 defines 
an area where the location of the complex poles of R. F. 

produce ay always in the second quadrant. This area is 

defined in terms of known parameters. It is that area which 


lies to the left of the complex pole of G,, where any complex 


1? 
root of R. F. will have a damping factor (& ) greater than or 
equal to that of G,, and will have a damped frequency of 
ose@illation Cw.) greater than that of the complex poles of 
G,. This area, in general, represents the possible locations 
of the complex pole of R. F. which will produce a faster 


response than G, but will not effect an increase in over- 


1 
shoot since the damping factor is at least as great as that 


for G,- Thus, the most probable location for placement of the 


Sc 








poles of R.F. lies in an area in which ¥Y is known to be a 


second quadrant angle. For complex poles of R.F. in this area, 
the determination of a positive ¥ will predict a faster 
response and a negative xy will predict a slower response 

of the desired function as compared to the representing 
function (R.F.). As the absolute value of y increases, the 
responses of the two functions become more dissimilar. 

Accuracy of the prediction of the desired response 

from the representing response iS gained with experience in 
using this technique. The steps utilizing the ''representing 
function technique" for determination of the conditional 
feedback compensation G, are summarized below. 

1. Specify the desired disturbance response of the 
compensated system with respect to the original 
basic system response. 

2. Utilize the root-locus method to determine a G, that 
will place the closed-loop poles of R.F. in a position 
such that desired specifications are approximately met 
and the R.FP. response iS approximately of second order. 

a. Place poles in that defined area in which 
Y 1S a second quadrant angle. 

b. Choose Wo of R.F. complex poles such that 
mi Cee) 


WG 1S equal to the desired improvement 
1 


in the response time. 
c. Pick a damping factor to meet desired overshoot 


requirement. 


oe 








a. 


4. 


6. 


Determine correlation parameters. 
: 4 r 
a. Compute Cy = W(R.F.) 
b. Compute ¥ using rule of equation (5-8). 
Make first prediction of desired function overshoot. 
a. Using second order approximation technique, 
estimate overshoot of R.F. 
b. Sum the estimated R.F. overshoot and Cy: 
Predict the desired function fall-time. 


a. Estimate fall-time of R.F. response as, 
| W 6G) 
fall-time (R.F.) = fall-time (G,) x W_CR.FS) 
(5-12) 

b. Predicted fall-time of the desired function 
response is the above estimate modified by r 
Positive ¥ predicts less fall-time and negativey 
predicts greater fall-time than that of the R.F. 
response. The degree of modification is propor- 
tional to the absolute value of v . 

Make second prediction of desired function overshoot. 

a. Increase first predicted value if ¥ is positive 
and decrease if negative. The degree of modifica- 
tion iS proportional to the absolute value of Tv. 

Compare predicted response with specified desired 
disturbance response. If satisfactory, proceed to 
verify response of the actual system including the 


G, conditional feedback compensator. 
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6. Determination of conditional feedback compensator. 
1 


1 set G,G, 


determine the specific G, compensation through root-locus 


The representing function, will be used to 


techniques in the manner outlined in Section 5. Observation 


of the response of the basic system to unit step disturbance, 


Fig. 6-1, shows it to have the following characteristics: 


Fall-—time: nOer) Sec. 
Overshoot: 10.3% 
Steady-state error: 0 


Number of oscillations: 1 
In order to effect a significant reduction in the track- 
ing error of the system due to transient disturbances, the 
transient characteristics of the representing function were 
arbitrarily chosen as follows: 
Pall-time: Decrease by 300% to 500%. 


Overshoot: No increase over that of 
basic system. 


Steady-state error: SMa . 


Oscillations: Essentially the same as 
for basic system. 


The procedure followed in using the root-locus method 
was to choose a G, which would effect an increase of Wi6G)), 
the natural frequency of the basic system, by an amount of 
300% to 500%. At the same time, the real roots of the basic 
system were compensated in such a way as to minimize the 
residues of the real roots of the resulting representing 


function (R.F.). Having done this, the R.F. approximated a 


second order system and its time response was generally known. 
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The simple correlations developed in Section 5 would then 
predict the desired function response. 
The basic system closed-loop transfer function (G,), was 


Paven in Section 3 as, 


37(s + 1.7)(€s + 20) 


©1 - (se s0d74)(s enema 14.53)(s + 22.63 = G14uSs) 


therefore, 


1 (s+1.74)(st+22.63+ j14.55)(s+22.63-j14.55) D, 


1+G,G, ~ Cst+1.74) (s*22. 63+ j14. 55) (s*22.63- 314.55) Dy + 


' * 37Cst1.7)Cst20) N. : 
(62a) 


where G., = dD. 


Several values of G., were utilized in an attempt to 
meet the time response requirements arbitrarily chosen for 
R.F. The following are representative of the compensations 
attempted. 

K 
aaa 


psn) —. . . » D> ApeenaemEnine tion, 


a lag function, 


K( stb) 


a , a> Dy a lecadweunic tion. 


KS . 2... ., a& d@Yivartivesnunctaen. 


The G, compensation which gave the desired position of 
the roots of the representing function (R.F.) was, 


we 21.06(€s+96)(s+200) 
2 (s+20)(st+1000) : 


a lag-lead compensating function. Fig. 6-2 shows the root- 
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imecns of the complex roome, of the R.F. utilizing thas G,. 


The complex roots are shown at 


s, = -80-j52 ; 5s 


i = -80+j52. 


2 


The real roots of R.F. were solved numerically and were found 


to be: 
S2 = = 115 (i Sea -20 ; Se = -1660 . 


The natural frequency of the basic system, Wi 6G,), is 26a? 
rad. per sec. The natural frequency for R.F., WACR-F.), is 
95.35 rad. per sec. This is an increase of 350% over Wi 6G,)- 


Section 5 has shown the form of the representing 


function to be Cy + Cem 


for the complex roots of R.F. determined above. The phase 


sin(52t +¥ ) 


angle (¥Y ) lies in the second quadrant since the complex poles 
are located in the area that insures a phase angle in the 
second quadrant. The steady-state error constant (C,) was 
approximated using equation (5-11) to be .079. Since the 
roots of this function are the same as those of the desired 
function, the latter function, to a second order approximation, 


-80t 


takes the form, C.e Sin€S2tat + tee 


5 
In order to apply equation (5-8) for determining the 
phase difference (%¥), the poles of Bae Gi and the complex 
poles of R.F. were plotted in Fig. 6-3 and ¥Y was measured 
aoetl 5 degrees. 
The R.F., to a second order approximation, should 
give a fall-time of about 350% faster than that of the basic 
system, or applying equation (5-12), 


Fali-time (R.F.) = 26.9 rad/sec x 0.0375 sec = .0O1 sec. 
: rad/sec 
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Since ¥ is positive and a relatively small angle, it.is 


expected that the fall-time of the desired function response 
will be slightly less than that estimated for the R.F., that 
meee slaghtly Wess than .01 sec. 

The overshoot of the R.F., to a second order approximation, 


was 0.008. This value added to the C, of 0.079 gives a first 


1 
prediction of desired function overshoot equal to 8.7%. 

Since ¥ was found to be positive, the overshoot of the desired 
response should be increased to some value slightly greater than 
that of the first prediction of desired function overshoot, that 
is, slightly greater than 8.7%. Thus, the predicted response 
of the desired function falls within the parameters selected 


for the new compensated system. 


The desired transfer function is, 


1 a 26 D.(D, : No? 
C1+F_)(1+G,G,) ~ (D_+N_) (DD, + N_D, + N_NZ) 
1+F 1+G,G, DotN, DoD. e N (Po + NGN> 
The roots, (DoD, SN ee NANo?> have been found to be, 


(s+1.712)(s+20)(s+1660) (s+80+ j52)(s+80- j52). 


Therefore, the desired transfer function is, 


s“(s+1.74)(s+10)(s+20)(s+1000)(s+22.63+ j14.55) 
S+1.74)(s+1.712) (st20)(8+22.63+ 314.55) (8422.63-)14,55 


(St22, 0 gama wo) 
Cs+80+ 352) Cs+80- j52)(s+1660) ° 


which after simplification is, 


x 


s“(s+10)(s+1000) 


Cst1.712)(st+1660)(S+80+ }52)(S#80-j5a) (6-2) 
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Using the graphical residue method on equation (6-2), 
the output due to a load disturbance applied as a unit step 
quG;, 


80t =1./12v | 


CG = .99@e0 Sin(52t.+ 25580) = 0009656 


(ween oo 


Fig. 6-1 shows a plot of this response. The predicted 
response characteristics together with those obtained from 


Pig. 6-1 are as follows: 


Predicted Ac tu au 
Fall-time: < .O1 sec .OO7 sec 
Overshoot: 8.7% + 8.6% 
Steady-state error: 0 QO 
Oscillations: -- 1 


The fall-time has been increased by 535% over that of 
the basic uncompensated system response. The overshoot has 
been decreased by 16.5% and no steady-state error results. 
This improvement in the regulation response did not result in 
any change in the desired original input-output response of 
the basic system. 

In order to check the accuracy of the correlation 
technique rules developed, the response of the representing 
function to a unit step input was determined using graphical 
residue methods, and was found to be, 


-80t 


C .0488 + .709e Ssin(a2t * 2: 306)ae 


Cae.) 


000745e7 1° /12t 4 yng Q71660t 
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ieeeeo-1 shows a plotvofithawe R. F. response. The difference 
in phase angles between the R. F. and the desired function 
memegual to (2.530-2.308) radians. Thus, the true value of r 
is equal to 0.222 radians or 12.7 degrees. The * obtained 
through the correlation rule, equation (5-8), was 13 degrees. 


The estimate of the R. F. C, was .079, while the actual value 


il 
of C, was computed to be .0488. The accuracy of the 
correlation rules used in the representing function technique 


has been shown to be more than sufficient. The time saved 


in uSing this technique for determining the desired conditional 


feedback compensation was considerable. 
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(eeeecoelection of diStwebance input. 

The external load torque on the antenna due to wind 
gusting is a function of time with three independent 
Variables, wind velocity, wind direction, and the tracking 
direction of the antenna. To evaluate the compensation, it 
is desirable to investigate its effectiveness in reducing 
the greatest expected disturbance error. Since the three 
variables affecting the disturbance output are independent, 
the maximum error will occur when the errors due to each 
Variable are maximum. 

To investigate the effectiveness of the compensation, 

a wind velocity profile having the greatest expected severity 
was used. From consideration of data pertaining to wind 

gust velocities“, a representative velocity profile was 
developed which would produce the maximum expected disturbance 
error. This velocity profile is Shown in Fig. 7-1. A steady 
25 miles per hour wind was assumed to exist prior to the 

gust. The wind velocity then increased smoothly to a peak 

of 60 miles per hour in one-half second and remained at this 
value for one-half second. The maximum gust velocity normally 
decreases at a slower rate than it increases during the gust 
build-up. Since only that part of the gust profile resulting 
in the largest disturbance error is of interest, the trailing 


edge of the actual gust profile need not be considered. 


a R. Gregg, Aeronautical Meteorology, pp. 269-274, 1930 
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the wind velocity input to the analo 
computer Simulation described in Appendix IV. 
Representative wind velocity profile. 


The wind velocity profile was develo 
straight line segments approximatin 


smoothly increasing wind 


utilize a 


@ 


Note 
Figure 7-1. 





The true wind direction was considered constant since 
wind direction changes are most likely to be negligible 
during the short interval of the gust being used. Since 
true wind has been considered constant during the one-half 
second of applied gusting, the change in relative wind 
direction is a function only of antenna rotation. The 
targets tracked by this antenna are earth-orbiting satellites 
which require relatively slow tracking rates. Thus, the 
change in relative wind due to antenna rotation is negli- 
gible during the one-half second gusting interval. During 
the gust, the antenna was considered to be in that position, 
relative to the true wind, which would cause the greatest 
disturbance torque for any given wind velocity. 

A plot of wind moment coefficient Core about the 
antenna gimbal center vs. angle of wind incidence has been 
prepared for the given antenna from wind tunnel tests made 
on a 1/30th scale model.” The data is reproduced in Fig. 7-2. 
The use of a ''spoiler,'' fixed to the back side of the antenna 
around the periphery of the dish, was used to reduce maximum 
torque loading on the antenna by approximately 50 percent. 

The moment about the antenna gimbal center (M_) was 


derived in Appendix III as 


MaezO1.9o C 2 
m 


V 
g Ww 


(lb-ft) . (7-1) 


Scuggenheim Aeronautical Laboratory, California 
Institute of Technology, Report 737, 1960. 
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where C= dimensionless moment coefficient about gimbal 
6 Center 


M8 = wWindewelloci tweety sec). 


W 


For a maximum velocity of 60 MPH and a maximum Ce 


of 0.189 (130° angle of wind incidence), the maximum Mo 


3 pat t. tines 


computed from equation (7-1) is 137.3x10 
moment, he: is the external torque oT! applied to the 


output azimuth drive shaft. 
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8. Evaluation of the conditional feedback compensation. 

The evaluation of the designed conditional feedback 
compensation is made on the basis of a comparison between 
the uncompensated disturbance-output response and the 
compensated system response. 

A step position disturbance input was used in Section 
6 to determine the design of the compensation since this 
input represented the most severe disturbance, though only 
theoretically possible. The compensation was designed to 
provide for an arbitrarily selected improvement in system 
response. These responses were Shown in Fig. 6-1. 

A comparison of system response to a step velocity 
disturbance input was made to closer approximate the system 
response to wind gusting. These responses are shown in Fig. 
8-1. 

In order to verify the degree of response improvement 
indicated as a result of the above two disturbance inputs, 
the system response to the wind gust profile of Fig. 7-1 
was obtained. The electronic analog computer simulation for 
this disturbance input was formulated in Appendix III. 
Appendix IV presents the analog simulation used to determine 
the system response to the wind gust. Both compensated and 
uncompensated system responses to the gust profile of Pig. 7-1 
are shown in Fig. 8-2a. 

In order to note the system response to a wind gust 
more severe than that derived in Section 7, and, in addition, 


one of less severity, wind gust profiles of one-half and 
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twice the time interval of the profile of Pig. 7-1 were 
applied. These responses are shown in Fig. 8-2b and Pig. 8-2c. 
The frequency response curves for the compensated system are 
Pyown in Pig. 8-3. 

Table I shows the disturbance response characteristics 
of the uncompensated and compensated systems for the 
disturbance inputs applied. The most significant parameter 
upon which to evaluate the compensated system is the maximum 
error caused by the wind gust. Table I shows this error, due 
to the wind gust of Pig. 7-1, to be 7.2 seconds of arc for 
the uncompensated system. With the addition of the conditional 
feedback compensation, this error was reduced to 0.028 seconds 
of arc. It is emphasized that this substantial reduction 
in antenna position error due to wind loading has been effected 


without change to the original input-output response. 
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9. Conclusions. 

This study has shown that there exists a compensation 
technique which can be successfully utilized to reduce 
antenna tracking errors due to wind disturbances without 
effecting a change in tracking input-output response. 

Section 2 introduced the conditional feedback compensa- 
tion method and has shown that this method permits design 
specifications on disturbance-output response to be met 
without change to the tracking input-output response. A 
specific conditional feedback compensation was developed to 
obtain a substantial degree of reduction of wind disturbance 
tracking error. The degree of disturbance error reduction 
is presented in Table I of Section 8. The representative 
wind velocity Eaof ile deveteume in Section 7 was applied as 
the disturbance input, and the resulting error decreased from 
7.2 seconds of arc without compensation to .028 seconds of 
arc with conditional feedback compensation. This represents 
a 257:1 reduction in disturbance error. In addition, this 
substantial error reduction was achieved with a faster response 
time. The response improvement of the compensated system is 
made more evident by noting that the other response parameters 
listed in Table I show substantial improvement over the uncom- 
pensated system. A significant disturbance error reduction 
might be expected when it is noted that bandwidth has been 
increased from 2 cps. to 98 cps. 

In order to determine the system response due to any 
given conditional feedback compensator, two root locus plots 


would have been required to place the roots of the closed loop 
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function. It was anticipated that several values of 
compensation would have to be investigated before a satis- 
factory compensator could be obtained. The representing 
function technique was originated for the purpose of expediting 
the above repetitive process. This technique was employed 

in the determination of the desired compensator. The represen- 
ting function technique proved to be accurate, expeditious, 

and simple to apply. A step-by-step procedure for its 


application is given at the conclusion of Section 5. 
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AEE ENDIA I 
COMPLEX ZERO COMPENSATOR 


An electronic compensator, Pig. I-1, has been developed, 
utilizing an operational amplifier, which will essentially 


cancel the complex poles of a cascaded transfer function. 





Figure I-1. Complex zero compensator schematic. 


iy - 1, + iy = () Ci =19) 
2 Ces (l-2y 
ee Z 
i, = — e Vy aa 
2 Ry 2 0 ag 
13 = C4V,s 
= Li 1 
— NV eS) 
1 Ry Ro 2 
7 eR> 
1 Ry + R, + R,R5C5s 
fee a 
2 -R, + R~ + R,R-CusS (1-3) 


1 2 Lee 


4G. A. Korn and T. M. Korn, Electronic Analog Computers, 
pp. 417-420, McGraw-Hill Co., 1956 
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2 a 
1. 72-— = 1 rs 
4 Re Ra 4 5 
i. = C2V,s 
E i i 
a- = ViC = + = + Cys) 
Ra 2 Re Ra 3 
ee 
e Re + Ry + R2R,C28 
: Ez 
Se (I-4) 
4 Re + Ra - R2R{C2s 
Substituting (1-2), (1-3), and (1-4) into (I-1) gives 
ec E 
Cles + a eee) 
I Ry + Ro + R,R5C5s Re + Ry + R2R,C35 de 


Letting R, = R, = R, = R, = 1 megohm, then equation (1-5) 


becomes 
e E 
C,es + SSS OC es FCO (I-69 
1 eee + C suite 2x00 4 peo 
6 2 12 
e (2x10 Cis + C,C.s x10 + 1) : _E 
_ 
2x10° + C.sx10'2 2x10° + C,sxl0-- 
2 : (ia 


If C, is made equal to C2, the denominators of (1-7) become 


equal, and (I-7) reduces to 


Ee. 2 2 z it 
eee 2° «SS 6+ GE C.c.xiole ? 
“ 2 
(I-8) 

If the capacitors are valued in micro-farads, equation 
(I-8) becomes 

E 2 Z if 

= = -C,C,(s" + st ) 

e LZ C, CC, (1-9) 
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Equation (I-9) is of the form of a constant times a complex 
zero, where C, and C, are picked so that the coefficients 
in this equation equal those of the complex poles it is to 
eancel. 

Depending on the values of capacitors used, it may be 
necessary to modify the circuit values slightly to suppress 
high frequency oscillations in the compensator amplifier. 
Reduction of the value of C, or placing a relatively small 
resistance (about 0.1 megolms) between C. and ground has 
been found effective. Such changes naturally reduce the 
cancellation effectiveness of the compensator, but to a 


degree generally within acceptable limits. 
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APPENDIX II 
GRAPHICAL RESIDUE TECHNIQUE FOR DETERMINING 
THE TRANSIENT RESPONSE EQUATION OR A SYSTEM 
Consider the disturbance-output transfer function 


of the uncompensated basic system: 


ewe i 7 s“(s + 10) - 
U~ T+ F ~ (s+ 1.74)(s * 22.63 - 314.55) 
x - (Tie 
$+ D2nCseeeeooy 
Let UCt) = A u(t), where A = 1. Therefore U(s) = = , aad 


equation (II-1) becomes 


s( SLO 


Suei(s + 1.74)Gs + 32.03 — Miss Geneon egI4s55) | 
(II-2) 
The poles and zeros of equation (II-2) are plotted in Fig. II-1l. 


The residue at (-1.74) is equal to 


TI'( vectors from the zeros to the pole ) e(pole)t 


TJ’ vectors from the remaining poles to the pole) 


or, (ry - z,)(r, - Z5) — 
[ry ao * jw) | Ir, = Ce jw) | (Tisep 


Prom Fig. II-1 equation (II-3) becomes 


(1.74 7 180° )(8226Gneae 1.74 /180° t 
a a a re e 
(25.45 /27 - 9,)(25.45 /@,) 

or -0.0222 e  - | aan (11-4) 
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The residue at the complex pole ( - 22.63 + j14.55) 


moeequal to: 


[«- X + jw) 2 Zz, | | (- a + jw.) - zp, | (GOs. jw. )t 


SS eee ae 
(- a + je - r, | Is oO + jw.) - (-@ - jw.) | 
| (Ill=5) 


From Fig. II-2, equation (II-4) becomes 
(26.9 age (19.25 /®>) 
(25.45/11 -8, (2 x 12 /@ 3) 


(26.9 ZY (19.25 [O>)  GPPne sme 1455) 
a 


2 (=22263 + 314255 


(25.45/T-0, (14.55) 25 
= 1.395 eme-- °° 53 ncia, 5a eee @, - IT+ Q,) 
1.395 eee °° 5 inCad Sone (Treen 


Combining equations (II-6) and (II-4) yields the complete 


transient response equation for the system Er Cie 
+ 


mee .O3T -1.74t 


weet 39D e sin(14.55t + 2.319) - 0.0222 e 


By this same technique, the other disturbance-output 
transfer functions used in this thesis can be converted to 
their equivalent transient response equations for various 
input Signals. A listing of these transient response 


equations for unit step and unit velocity inputs follows. 
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il 5“(s + 10) 
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i+ F_ PA ” Cs + 1, 7) Coe 63 — Sra. oer aie. SS) 
Ct) = 1,305 6 22: 93t gaia .55t + omamlo) - 0.0p00 © 1 ae 
oan = 0.0520 e ~°°°9°! Gantiayss5t — ©2247) + lonolle7 « ous 
E 1 
U" I+ GG, 
—=(s + 1.74)€s + 22.63 - j14.55)(s * 22.63 + j14.55)(s + OG 
Cs + 1.412)Cs + 80 - j52)Cs + 80 + j52)Cs + 1660) 
cit) = 0.0488 + 0.709 e ~849* sincsat + 2.308) + .000745 e ~t-/14t 
+ 0.426 e ~1600t 
8 = 0.00172 + 0.0488t + G100728NeE UlenCSet =o loaEn 
+ 0.000423 e ~2°/1¢* _ 9, 900257 e ~1060t 


—_-_ -_-__-]- S>-llCelleESll Slee me mm me i mm mm meee ieee Peete 


s“(s + 10)(s + 1000) 


Comm ievi2)(s + 80 = j52)Cs + 80 + j52)Cs + 1660) 
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AP eeNDLTA eit 
DEVELOPMENT OF DISTURBANCE TRANSFER FUNCTION 


The purpose of this Appendix is to develop the relationship 
between the disturbance input (U), a function of the disturbing 
wind velocity CF) and the antenna position angle (0, ). The 
development of this transfer function was undertaken in two 
steps. First, the transfer equation of the disturbance torque 
(TD with respect to Le was determined. The transfer function 
of T, to the antenna position angle (O, ) was then determined 
and the two steps combined to give the transfer function of ur 
to QO; - 

Fig. 7-2 presents a plot of the non-dimensional moment 
coefficient (Coe? of the antenna as a function of relative 


wind direction. Using Fig. 7-2 to determine Cag? the moment 


about the antenna gimble center (M_) which is also the external 


torque (T.), was found as follows:? 
Toe = Cn Se (III-1) 
_ a 2 . 1b 
where, q=28V_, the dynamic pressure (zF2), 


5 = 2630 ft~, antenna dish frontal area, 


§ = air density (Abesec*), 


yr = 30 ft, antenna dish radius. 
The data of Fig. 7-2 was taken at a temperature of 20°. 
and a barometric pressure of 742 mm Hg. 


. * . a = ee, 
The air density (f ) is defined as, f ERT ETI => 


SJ. H. Dwinell, Principles of merodynamies, pp. 25-26% 
51-52, McGraw-Hill Co., 1949 
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where, p = absolute pressure ( 5 
or 
= fit ; 
g =z 2 5 eecelération due to cravity 
sec 
R = 25.3 na , gas constant 
OR 


T = absolute temperature (°R) . 
The temperature of 28°C is equal to an absolute temperature 
of 541.9°R. A pressure of 742 mm Hg. is equal to 2075 ibs/ft 
at 289°C. Substituting these constants into equation (III-2) 


gives, 
P2902 235 ib-See>/ tee 


Solving femeqg in terms of Ve and substituting into 
equation (III-1) gives, 


¢ eO1.9 C_. Vo eee (ITT 
€ mg w 


Equation (3-18) of Section 3 shows the relationship 
between the error signal (X) and the external load torque 
ar with respect to the position output angle (0, ) of the 
hydraulic motor. Fig. 3-2 shows the hydraulic system configu- 
ration. Thus, the actual input to the hydraulic function (G,) 
of Pig. 3-2 for an applied external wind torque is the sum 
of the error signal (X) and some function of external torque. 


The hydraulic motor output angle (@,) can then be written as, 


0, = G, xX + GG, fT, CIT tay 


where G, is some transfer function which carries external 
torque into an error position to be applied to the hydraulic 
transfer function (G,)- Letting X in equation (III-4) be ecauman 
to zero, then e, = G,G, - ; CI Tiley 


TD 





Substituting X equal to zero into equation (3-18) gives, 


O, (-) K, (L+gs) 
cee dn NJ, 8°+£, stK,) [dm s+(L+ps)J_s“+(L+gs)K, (J, s°+f,s) 
(Tilia) 


Substituting equation (III-6) into (III-5) gives, 


alen 


V 
(-) Ky (L+Es) 





2 a V 2 vo VK 
S | am NCJ, s +f, stK, )+(Ltgs)NJ_s(J,s *£ Stk, )+ (Legs) LCI, sf, ) | 


(III-7) 
Making the same aSsumptions as in Section 3, that is, Ky is 
fL 
fey large and@, ——__ = 0), equation (1@l-7) ‘simplit lec geen 
ec) +J.) 
m  & 
Ne 


(ae) (L+%s)NB/V(N“J +J_) 
= B Me 
G,G, = Sas) or (Ilil=3e5 
S s +BLst+tB dmN 
Veowvenert at) 
VC(N Jat 3p 
Dividing equation (III-8) by G, as Piven equation (3-21) 


gives, Vv 
(G,G,) -(L+gs) 


G2 ee (III-9) 
h Vv 


Fig. III-1 shows the final block diagram of the hydraulic 
system with the added complex zero compensation (G) and 
tachometer feedback (G,-) described in Section 4. In addition, 
the application of a through the synthetic transfer function 


G, is represented. 


il 
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Figure III-1. Block diagram of final hydraulic system. 


Assuming that @, is due to the external torque alone, 


iL 
ore. X' = 0), then 


= = “1h (III-10) 
T N(1 + G;-GG,) ’ 7 
e £269 
where, Gr = Kes, the tachometer feedback 
Go = Ko x (complex wcots on Gi); the complex zero 


compensator added in Section 4. 


Substituting equation (III-8) into (III-10) and simplifying 


gives, 
Or (s+BL/V) / (NOS +J,) 
Ne ; S s°+BL/Vst+ B dm2 Ne 1s K¢KoKyBN~ dm 
V(NJ_+J,) VJ, +N°S_) 


(II l= 
Substituting the hydraulic flow coefficients listed in 
Section 3 and the gain constants developed in Section 4 into 


equation (III-11) gives, 


L  £.5ixl0 Come rad/lb-in. (1II-12) 
ee 
T. s(6°+4264 3008 


oT 





Combining equations (III-12) and (III-3) enables the 
solution of the output position angle of the antenna for 


any given wind velocity. Thus, for a constant .. 


g 9 
Q -§8 
L - 201.9 Cng Klee x 1 0 (st+42) (III 13) 
Z 2 . s 
Vw s(s~ + 42s + 3228) 
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APPENDIX IV 


ANALOG SIMULATION OF CONDITIONAL 
le EPBACK “COMPENSATED SYSTEM 


Fig. 5-2 shows the block diagram of the complete 
compensated system under study. It has been shown that only 
the response of the system to a disturbance input is of 
interest since the response to the tracking input (R) is 
known to remain unchanged with the addition of the conditional 
feedback compensation. Fig. IV-1 shows the block diagram of 


the compensated system with the tracking input equal to zero. 





Figure IV-1. System with conditional feedback and R = OQ. 


The transfer blocks of Fig. IV-1 were defined as follows 


m7 section 4: 
mes + 1.7) (Sez 
5° (s + 10) 


21.06(s + 96)(s + 200) 
a (s + 20)€s + 1000) ’ 


i 
O 


G 


and where U is the disturbance input, (O,), the antenna 
position due to an applied wind velocity. Equation (III-13) 


has shown the transfer function from wind velocity (Vv) to 


fee 





=. ombe , 


U_ _ OL _ 201.9 Cmg x 1.51 x 107%(s + 42) 
View Vay ee + 42s + 3228) 





The analog simulation of the compensated system of 
Fig. IV-1 was undertaken in order to determine the effects 
of the added conditional feedback compensation on the system 
response to a simulated wind velocity disturbance profile. 
The following is a list of equipment utilized in the 
Simulation of the compensated system: 


3- Donner analog computers, Mod. 3000, 10 channels 
Dem COMpUutEer. 


1- Brush dual-channel D.C. amplifier. 
1- Brush dual-channel pen recorder, Mod. BL-202. 


1- Donner variable-base function generator, Mod. 
3750 


1- Donner electronic function multiplier, Mod. 3730. 
1- Donner null. voltmeter, Mod. 5002. 


1- Hewlett Packard low frequency function generator, 
Mod. 202A. 


1- Moseley Autograph x-y recorder, Mod. 2-S. 
The simulation of the disturbance input function was 
undertaken in two parts. First, the disturbance torque (TD) 


was simulated using equation (III-3), 


T = 201.9C Vo" Sens 
e mg Ww 


The representative wind velocity profile of Fig. 7-1 
was Simulated as a function of time using a variable-base 


function generator. The simulated profile was generated using 
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a step velocity input as a time base and an eleven straight- 
line segment approximation of the wind velocity profile. The 
constant initial velocity of 25 MPH was then added to the 
output of the function generator. The resultant was the 
Simulated profile of Fig. 7-1. This velocity profile was 
Squared using the function multiplier. Sanece Cag was considered 
constant during the interval of the gust, the squared velocity 
was multiplied by the appropriate constant to give 1 & 
Equation (III-12) gives the transfer function between 
T, and the output position due to wind loading (QO, ). The 
analog simulation of this transfer function requires an 
integrating output stage. This equation represents an open- 
loop transfer function. Therefore, there is no feedback loop 
eround this owfput stage. Since the input to this inteerarer 
is non-zero, this integrating stage would saturate. The 
Simulation was modified to include this integrating stage 
within the feedback loop of the servo system, resulting 
in an output from the disturbance simulator of the velocity 


(QO, ). Modifying equation (III-12) gives, 


BL _ K Cs 7 eee ) 
T — 
S (5° 


Un 
= = 


e + 425 + 3228) 


where Ie 1S equal to 1.5islon 4 


The velocity 9, was added to the velocity output of eae 


I 
transfer function F and integrated to obtain the system 
Posit@on OUtpine)(C). 


Pig. IV-2 shows the complete schematic of the analog 
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computer simulation used. Values of components used as well 


as the primary outputs are indicated on the figure. 
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